Energy loss of fast quarks in nuclei 
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We report an analysis of the nuclear dependence of the yield of Drell-Yan dimuons from the 800 
GeV/c proton bombardment of 2 H, C, Ca, Fe, and W targets. Employing a new formulation of the 
Drell-Yan process in the rest frame of the nucleus, this analysis examines the effect of initial-state 
energy loss and shadowing on the nuclear-dependence ratios versus the incident proton's momentum 
fraction and dimuon effective mass. The resulting energy loss per unit path length is —dE/dz = 
2.32 ± 0.52 ± 0.5 GeV/fm. This is the first observation of a nonzero energy loss of partons traveling 
in nuclear environment. 
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For many years it has been suggested that quark en- 
ergy loss might give rise to a nuclear dependence of 
the cross section of Drell-Yan || (DY) production. When 
a proton enters a nucleus the first (soft) inelastic collision 
liberates a quark, which then loses energy via hadroniza- 
tion (due to confinement) and interaction in the nuclear 
medium. A lepton pair created from a subsequent inter- 
action then has reduced energy compared with the DY 
process on a free nucleon. The goal of the present anal- 
ysis is to search for this effect in the nuclear dependence 
of the DY process. 

Fermilab E772 made a precise measurement of the nu- 
clear dependence of the DY process using 800 GeV/c 
protons. The experimental details of E772 have been 
described previously |^-^|. Briefly we indicate those 
germane to the present discussion. Muon pairs were 
recorded from targets of 2 H, C, Ca, Fe, and W, in the 
mass range M > 4 GeV/c 2 . Excluding the T resonance 
region, 9 < M < 11 GeV/c 2 , we reconstruct 2.5 x 10 5 
DY dimuons. The spectrometer acceptance for this sub- 
set of the data had transverse momentum coverage out 
to 3.5 GeV/c. Since E772 was designed to make a preci- 
sion comparison of the yields of dimuons from the heavy 
targets to that from 2 H, relative target-to-target normal- 
ization errors were kept to < 2%. 

The parton model description of high-energy processes 
is reference- frame dependent. Because energy loss is most 



commonly described in the rest frame of the nucleus, it is 
best to adopt this frame for the description of the DY pro- 
cess as well. In the target rest frame the DY process for 
proton- nucleon collisions is treated as bremsstrahlung || : 
An incident quark with momentum fraction x q emits a 
virtual photon that carries a fraction x\ = x\/x q of the 
quark momentum. The inclusive cross section for the 
production of lepton pairs with momentum fraction x\ is 
given by 



dxi 



dXqFq{Xq) ~ 



dx\ 



(1) 



where F^(x q ) is the quark distribution function of the 

proton and da'^(x1, M 2 )/dx\ is the quark-nucleon dif- 
ferential cross section for lepton-pair production [j^— [lT|] . 

Nuclear effects modify this in two important ways. The 
first is the possibility of quark energy loss in the nuclear 
medium - the main subject of this manuscript. The sec- 
ond is shadowing, a phenomenon well known from nu- 
clear dependence studies Q of a closely related process, 
deeply- inelastic lepton scattering (DIS). Energy loss and 
shadowing are shown pictorially in Fig. ^. Since the two 
processes produce apparently similar effects in proton- 
nucleus collisions, it is necessary to adopt a consistent 
analysis where both are considered on the same footing. 
The framework for accomplishing this is detailed in the 
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following paragraphs, first for energy loss, then for shad- 
owing. 

Consider a proton entering a nucleus (Fig. Q). The 
first inelastic interaction, at point z\ , removes the coher- 
ence among the soft projectile partons, which then move 
apart losing energy as they would in the vacuum. A 
quark continues to propagate arriving at point z^ where 
a DY interaction takes place with diminished energy 
i q E p — x q E p — AE, where AE is the energy loss to 
be measured and E p is the energy of the proton beam. 
Correspondingly, one has x\ = x\/(x q — AE/E p ). As- 
suming the rate of energy loss is constant, AE oc L, 
where L = zi — z\. Due to energy loss, the ratio of p-A 
to p-N cross sections versus x\ is 



Ixi+AE/Ep Fqi x q) 



qN 
1 DY 



{x\,M ) 



(2) 



A distribution of energy losses occurs. We have cal- 
culated this distribution using Glauber theory. In the 
calculation, the relatively small probability for inelastic 
collisions of the incident proton leads to a significant re- 
duction of (L) = (z2 — Z\) with respect to the mean path 
length Lq. For example, for tungsten we find (L) = 2.4 
fm, whereas for a uniform sphere Lq — SRqA 1 ! 3, / A = 4.9 
fm. 

Shadowing is the well-known reduction of the cross sec- 
tion per nucleon, observed experimentally in DIS |12| for 
x less than about 0.07. In proton-nucleus DY production 
a reduction in the cross section per nucleon at small x is 
seen at the highest available proton energy, 800 GeV/c 
P, fb3|Jil|l . It is, however, an open question whether this is 
shadowing, energy loss, or both. The framework for an- 
alyzing shadowing in the rest frame of the nucleus, given 
below, is crucial to resolving this puzzle [p|, pr|jr^Jl^1 . 

In both DIS and DY shadowing occurs when the nu- 
clear coherence length grows larger than the distance be- 
tween nucleons, « 2 fm. The coherence length is a mea- 
sure of the lifetime of the fluctuation of a quark into a 
virtual photon and residual quark. For the DY process, 
the mean nuclear coherence length is given by 



lr 



2E q x q 1 {\-x\) 



(1 - x\) M 2 + (x\ m q ) 2 + k 2 



(3) 



where E q — x q E p and m q are the energy and mass of the 
projectile quark which radiates the virtual photon. The 
resulting lepton pair has an effective mass M, a trans- 
verse momentum k?, and carries a fraction x\ of the ini- 
tial momentum of the quark. The mean coherence length 
for the kinematic conditions of E772 has been evaluated 
in Ref. p"7j ] by integrating over x\ and kj*. This is sim- 
ilar to a previous very successful treatment of shadow- 
ing in DIS pq ]. The result is shown versus x\ in Fig. ^ 
for various fixed values of x-±. The coherence length is 
nearly independent of M 2 at fixed xn , but it vanishes at 
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FIG. 1. Schematic representation, in the nuclear rest 
frame, of two processes producing a nuclear dependence of 
the DY cross section, (a): Energy loss. A proton entering 
a nucleus undergoes its first inelastic collision at point zi, 
liberating a quark. The quark propagates, losing energy to 
hadronization (dashed lines), to point Z2 where it undergoes a 
hard (DY) interaction, producing a virtual photon that decays 
into a lepton pair, (b): Shadowing. A fast quark undergoes a 
virtual fluctuation into a photon and quark. The quark prop- 
agates into the nucleus, liberating the fluctuation. Since the 
coherence length is large, the entire nucleus participates as a 
single entity. 



xi — > 1, violating factorization. We note that the val- 
ues of l c given by Eq. || are significantly smaller than the 
commonly cited l c « \j2m^X2 (see discussion in fl6f|). 
In the weak shadowing approximation [lql. 



R'Jffi{x lt M') 
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Veff {T A )F\(q c ) 



(4) 



This expression is accurate for the whole kinematic range 
of E772. Here, (Ta) is the mean value of the nuclear 
thickness function, q c = l/l c , and 



Fl(q cl 



1 



A(T A ) 



d 2 b 



dz i 



' PaQ>,z) 



(5) 



is the longitudinal nuclear formfactor where PaQ>, z) 
is the nuclear density. The effective cross section is de- 
fined as H a e ff — (a 2 {x\tt)) I Mx\tt)) , where g(x\tt) 
is the qq dipole cross section |20(], and tt is transverse 
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FIG. 2. The mean coherence length as function of x\ at 
fixed values of xi — 0.02, 0.03, ... 0.08, evaluated for the 
kinematic conditions of E772. 



separation between the virtual photon and the quark 
(Fig. |); a eff is in the range 3.5-5.5 mb for E772 @. 

Energy loss and shadowing provide mechanisms for nu- 
clear suppression in the DY process that are effective 
in different regimes. If l c is short, no shadowing oc- 
curs (F\(q c ) — > in Eq. ||), but energy loss can reduce 
the yield of DY pairs. In the opposite limit, l c 3> R A , 
shadowing achieves its full strength (F\(q c ) — > 1). Here, 
initial state interactions do not affect the DY cross sec- 
tion except for transverse momentum broadening p],pl|. 
Our ansatz is that the transition between these limit- 
ing regimes is controlled by F\(q c ). The only expres- 
sion which is linear in F A (q c ) and has the right limits at 



-» and q c — > oo reads 
R A/N (x u M 2 ) = (i^^O-l) 



(6) 



Given the above framework, the energy-loss analysis 
was accomplished as follows. The quark- nucleon cross 
section in Eqs. |l| and || was assumed to have the form 



daf Y (M 2 ) 



dx\ 



K(M 2 ) x (1 - x\) 



q\m 



(7) 



where K and m were determined from a fit to the p± 2 ff 
data with Eq. 0. It was found that m did not change 
significantly with mass bin; thus one slope parameter 
was sufficient to characterize the full range of the p+ 2 H 
data. Because the mass-dependent normalization (K) of 
the q — N cross section occurs in both numerator and de- 
nominator of Eq. ^, the energy loss term in Eq. |^ becomes 
independent of mass. 

Thus, in this formulation, energy loss and shadowing 
have different kinematic dependence. Therefore is essen- 
tial to analyze nuclear-dependence ratios that are binned 
in dilepton mass. This allows l c to be evaluated for each 



bin in mass and x\. A two-parameter fit using Eq. || 
was applied to the heavy-target cross section ratios for 
C, Ca, Fe, and W; the parameters were —dE/dz and 
an overall normalization factor, C. The systematic nor- 
malization error, ±1%, was treated as an additional sta- 
tistical error. The fit yields a substantial energy loss, 
-dE/dz = 2.32 ±0.52 ±0.5 GeV/fm (statistical and sys- 
tematic), with C = 1.010 ± 0.006, consistent with the 
E772 normalization uncertainty The systematic er- 
ror associated with —dE/dz arises from uncertainties in 
cut parameters, the range of applicability of Eq. [7], and 
the shadowing analysis (discussed below). Fits to W/ 2 H 
and to C / 2 H in four mass intervals are shown by solid 
curves m Figs. | and |. 
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FIG. 3. Ratio of tungsten to deuterium Drell-Yan yields 
per nucleon versus X\ for different intervals of M. Dashed 
curves correspond to net shadowing, solid curves show the 
full effect including shadowing and energy loss. 

Our the analysis depends critically on having separated 
the effects for energy loss and shadowing. Dashed curves 
show the net shadowing contribution. Nuclear suppres- 
sion of the DY cross section for tungsten is mainly due 
to energy loss. On the other hand, energy loss effects for 
carbon are small, the main contribution to nuclear sup- 
pression arising from shadowing. This difference between 
the A- and M-dependence of energy loss and shadowing 
permits the two effects to be disentangled. 

We have checked the sensitivity of the results by the 
following tests: (i) eliminated shadowing in Eq. || by 
fixing F\ = (-dE/dz = 2.24 ± 0.53); (ii) doubled 
the shadowing corrections, 1 - R s A , a ^ =>■ 2(1 - R s A h / a ^) 
(—dE/dz = 2.64±0.53); (iii) mixed energy loss and shad- 

n A/D X 



owing effects differently from Eq. ^, R A /d 
R A/ a D (-dE/dz = 2.35 ± 0.53: This is the E866 pro- 
cedure [0.); (iv) selected for analysis only data with 
small 22 < x™ ax : Within error bars —dE/dz is constant 
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FIG 
Fig. |. 
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4. Ratio of Carbon to deuterium. Same labeling as 



for 0.3 > xf ax > 0.12. Thus -dE/dz is subject to only 
small variation within the error bars. These modifica- 
tions all lead to a significant growth of x 2 ■ 

Recently the E866 |L3) collaboration analyzed DY 
nuclear dependence data from targets of Be, Fe, and 
Cu. The E866 data set, 1.3 x 10 5 muon pairs in the 
range 4 < M < 8.5 GeV/c 2 , was more concentrated 
at low xi than the present one. E866 subtracted a 
phenomenological shadowing contribution, yielding new 
"shadowing-corrected" nuclear dependence ratios. Shad- 
owing was calculated employing the results of the global 
phenomenological analysis of DIS and DY data by Eskola 
et al. |22j (EKS). The EKS analysis itself included the 
DY data from E772, with the presumption that the low- 
X2 nuclear dependence arose entirely from shadowing. 
This clearly introduced an inconsistency into the E866 
search for energy loss. Considering the critical impor- 
tance of separating shadowing and energy loss, which, at 
800 GeV, can be achieved only via mass-binned nuclear- 
dependence ratios, it is not surprising that the E866 anal- 
ysis missed the effect. 

The value of dE / dz determined here is not very differ- 
ent from that found many years ago by Gavin and Mi- 
lana (GM) Q| using the mass-averaged W-D ratio from 
E772. In our model shadowing is a very small effect for 
the mass-averaged W data, so the GM analysis, which 
ignores shadowing, should not be too far off. The GM 
value, dE/dz k, 1.5 GeV/fm, should be increased by a 
factor of « 2 to account for the reduction of the effec- 
tive nuclear path length, discussed earlier in connection 
with Eq. (||). Unlike the GM model, our analysis pre- 
sumes a constant dE/dz, yielding an energy loss that is 



independent of laboratory beam energy (see Ref . [£3| ) . 

Much theoretical attention has been devoted in re- 
cent years to the elucidation of the QCD analogue of 
the famous Landau-Pomeranchuk-Migdal [^4| effect (see 
Ref. (^5|). Gluon radiation induced when a quark pene- 
trates nuclear matter leads to additional energy loss pro- 
portional to the square of the path length traversed. Us- 
ing the measured transverse-momentum broadening Q| 
this can be estimated for Tungsten as rising to a maxi- 
mum value, — (dE / dz) ra d ~ 0.2 GeV/fm. Thus for cold 
matter radiative energy loss is not a large contribution 
to the total. 

In light of the present finding that quark energy loss 
is significant, one should re-examine the role of energy 
loss in the nuclear dependence of J/tjj production (see 
Ref. (^6|). It was demonstrated many years f|,|27j ago 
that the ^-dependence of J/ip suppression at energies 
150 — 300 GeV could be well described by energy loss, 
with —dE/dz pa 3-4 GeV/fm. The larger value for a 
projectile gluon is consistent with an enhancement due 
to the Casimir factor, 9/4. 

In summary, we have made the first determination 
of quark energy loss using an analysis that takes into 
account nuclear shadowing. The result, —dE/dz = 
2.32 ± 0.52 ± 0.5 GeV/fm, is in approximate accord with 
the theoretical expectation that energy loss should be at 
least the order of the string tension, k w 1 GeV/fm. At 
800 GeV mass binning of the nuclear dependence ratios 
is crucial to the separation of energy loss and shadowing 
effects. It would be very desirable to have precise mea- 
surements of the nuclear dependence of DY production 
at lower beam energies (100-300 GeV) where shadowing 
disappears and energy loss would provide the dominant 
nuclear dependence. 
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